Observation of attraction between dark solitons 
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We demonstrate a dramatic change in the interaction forces between dark sohtons in nonlocal 
nonlinear media. We present, what we believe is the first experimental evidence of attraction of 
dark solitons. Our results indicate that attraction should be observable in other nonlocal systems, 
such as Bose-Einstein condensates with repulsive long-range interparticle interaction. 
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It is commonly accepted that solitons, i.e. localized 
waves propagating without changing their shape, are 
ubiquitous in nature 1] and are native to many diverse 
systems like plasmas, molecular chains, spin waves, at- 
mospheric physics, superfluidity, nonlinear optics, and 
Bose-Einstein condensates (BECs). A soliton can form 
when the dispersive or /and diffraction processes associ- 
ated with the finite size of the wave are counterbalanced 
by the wave self-induced change of the properties of the 
medium. In the context of BECs, for example, the soli- 
ton represents a coherent excitation of a matter wave Ij , 
whereas in optics it is a localized light beam or pulse . 
The continuous interest in solitons is stimulated by their 
unique collisional properties, i.e. they behave like par- 
ticles displaying forces in their mutual interaction. Fur- 
thermore, the fundamental features of their interaction 
are of rather universal character. Thus, for example, 
matter wave solitons interact basically in the same way 
as optical or plasma solitons. 

There are two fundamental types of solitons: bright, 
in the form of a localized structure ^ and dark, in 
the form of a localized dip on a plane-wave back- 
ground H, IE S S ^ • Already early studies pointed 
out a distinctive difference between the interaction of 
bright and dark solitons. Bright solitons can either at- 
tract or repel depending on their relative phase. The 
phase between dark solitons is fixed as they are formed 
on a single background wave and they are believed to al- 
ways repel [3, 0, 0, This fact imposes a funda- 
mental limit to the applicability of dark solitons and all 
existing soliton applications are currently based only on 
bright solitons. 

Recently, we predicted theoretically |3 that the na- 
ture of dark soliton interaction may change drastically in 
the presence of a nonlocal response of the material. In 
this letter we report the first experimental observation 
of attraction between spatial dark solitons in any physical 
system. Our findings open new possibilities for control of 
the interaction between dark solitons, which we believe 
will revive the interest towards them and will allow for 
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FIG. 1: (color online) (a) Trajectories of two attracting dark 
solitons separated at GO/xm. (b) Initial intensity (solid line) 
and corresponding index structure (dashed line) generated 
by the nonlocal nonlinear response, giving rise to attractive 
(dashed arrows) and repulsive forces (solid arrows) on the 
solitons. (c) Trajectories of two closely spaced dark solitons 
(20/im) with dominating repulsive interactions. 



their broader applicability. 

In nonlocal nonlinear media the nonlinear response in- 
duced at a certain point is carried away to the surround- 
ing regions. In this way a narrow localized wave can in- 
duce a spatially broad response of the medium . Spa- 
tial nonlocality is an inherent property of many physical 
systems. It often results from transport processes, such 
as atomic diffusion 16] or heat conduction 0. Spatial 
nonlocality is also natural for media with a long-range 
inter-particle interaction including, for instance, dipolar 
BECs 18j, or nematic liquid crystals with long-range 
molecular reorientational interactions ^2|- It appears 
that nonlocality leads to novel phenomena of generic na- 
ture. For instance, it may promote modulational insta- 
bility in self-defocusing media [13,1^, as well as suppress 
wave collapse and stabilize multidimensional solitons in 
self- focusing media [23, 1^ • Nonlocal nonlinear- 
ity may even describe parametric wave mixing 26] . Fur- 
thermore, nonlocality significantly changes bright soliton 
interaction 27]. 

To test experimentally how the nonlocality affects the 
forces acting between dark solitons, we considered the 
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propagation of an optical beam in a weakly absorbing 
liquid. Light absorption increases the temperature of the 
liquid and subsequently decreases its density and refrac- 
tive index, resulting in a defocusing nonlinearity. In ad- 
dition, heat conduction leads to a temperature and con- 
sequently refractive index profile much wider than the 
light beam itself, indicating the inherently nonlocal char- 
acter of the thermal nonlinearity. The choice of the liquid 
medium provides additionally the advantage to monitor 
the beam profile along the entire propagation without 
destructive intervention to the material. 

The propagation of an optical beam along the 2;-axis 
of a weakly absorbing liquid is described by the system 
of normalized equations [13 for the slowly varying field 
amplitude 2:, t) 
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Here the longitudinal z and transverse x coordinates are 
scaled by the wave vector k = 27rno/A {z kz/2^x 
kx). The parameter a (scaled as a a/k) represents 
the linear absorption. Time in Eq. is normalized to 
k^D^ where D is the thermal diffusivity of the medium, 
and A is the laser wavelength. Eq. (pP) is a heat equa- 
tion describing the temporal and spatial dynamics of 
the refractive index change of the medium An{I) in- 
duced by a heat source in the form of a beam with 
intensity I{x^z^t) = 2:, Eq. (|2j is a nonlin- 

ear Schrodinger equation governing the evolution of the 
beam amplitude in the presence of nonlinear refractive 
index change. We solved these equations numerically us- 
ing a split-step Fourier method for the propagation equa- 
tion and a finite difference method for the heat equation. 
For simplicity we assumed unrestricted heat flow in the 
transverse direction. Using typical thermal parameters of 
mineral oils yields the diffusion coefficient D = l{)~^m? / s 
and linear index no = 1.5 for A = 0.532/im. The absorp- 
tion coefficient is set to a = O.Olcm"^ which results in 
10% power loss over 50mm of propagation. As an ini- 
tial condition we used a broad Gaussian beam of peak 
intensity /q and full width at half maximum (FWHM) 
2.8 mm, with two closely placed tt phase jumps. Such 
initial conditions result in the formation of two "black" 
solitons with a zero individual transverse velocities E^. 

An example of steady- state trajectories of two adja- 
cent solitons separated by a distance of 60/im is shown in 
Fig.^a). The inward bending of the trajectories clearly 
indicates the presence of attractive forces. For longer 
propagation the separation between the two solitons os- 
cillates, leading to formation of an oscillatory-type bound 
state. This behaviour strongly differs from the interac- 
tion of dark solitons in local nonlinear media, where the 
soliton trajectories diverge due to repulsion. 
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FIG. 2: (color online) (a) Calculated output separation 
(z=50 mm) versus initial separation between two dark soli- 
tons for low (/o = 0.01 - squares) and high (/o = 0.5 - circles) 
intensity of the background beam. The different regimes of 
interaction are marked. (b,c) Output intensity (bottom) and 
index (top) profiles for two different initial separations - (b) 
repelling solitons (20 /xm) and (c) attracting solitons (90 /xm). 



The physics of this interaction can be intuitively ex- 
plained as follows. The two close dark solitons with an 
intensity profile depicted in Fig. Q^b) always try to repel 
because of local refractive index drop in the overlapping 
region (repulsive forces are indicated as solid arrows). 
However, in a nonlocal medium these solitons induce also 
a large scale change in the refractive index in the form 
of a broad trapping potential [Fig. Q^b, dashed curve)]. 
This potential provides an attractive force (indicated by 
dashed arrows) which counteract the natural repulsion 
of the solitons. Ultimately, the interplay between these 
two forces will determine the outcome of the interaction. 
While repulsion prevails for close initial separations and 
the solitons diverge as depicted in Fig.^c), its strength 
decreases with increase of the separation, allowing the 
nonlocality mediated attraction to become dominant at 
larger separations. 

The spacing between the solitons obviously determines 
the shape of the resulting potential, and subsequently the 
strength of their interaction. In Fig. EJa) we present the 
numerically determined separation between the two dark 
solitons, after propagation of 50 mm, as a function of 
their input spacing for low (linear regime - squares) and 
high intensity (nonlinear regime - circles) of the back- 
ground beam. For initial spacings larger than 180 jiia the 
final separation is close to that of the linearly diffracting 
dark beams as the nonlocality does not contribute appre- 
ciably to the already weak interaction. The most inter- 
esting regions in Fig. Elja) correspond to closely spaced 
solitons (20-180 jim). It is evident that in this region 
the separation between the solitons is a non-monotonic 
function of the initial spacing, being smaller than that 
in the linear regime. This behavior is a direct manifes- 
tation of the nonlocality induced attractive forces. For 
closely spaced initial dark notches (20-50/im), the gener- 
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FIG. 3: (color online) Experimental setup: M - mirrors, GS 
- glass slides, BE - beam expander, T - imaging telescope, P - 
pellicle, IL - imaging lenses, A - attenuator, NLM - nonlinear 
medium, and CCD 



cameras. 



ated solitons strongly repel and their final separation is 
larger than the initial one. In the case of intermediate 
separations (50-180/im), the attractive forces can balance 
the natural repulsion and the solitons exhibit oscillating 
trajectories. Plots in Fig.I2fb,c) illustrate the output in- 
tensity and the refractive index profiles (which plays the 
role of confining potential) in those two cases. For the di- 
verging solitons the potential has a form of two distinct 
wells separated by a barrier, whereas attraction occurs 
when the potential represents a single well. 

To investigate experimentally the interaction of dark 
solitons in a nonlocal nonlinear medium we used the ex- 
perimental setup shown in Fig. |31 A laser beam from a 
frequency doubled solid-state laser (Verdi-V, A = 532nm) 
was expanded by a system of lenses and passed through 
two closely overlapping microscope glass slides. The 
slides were subsequently imaged by a telescope onto the 
input face of a 50 mm long glass cell. The resulting beam 
diameter at the cell's input was FWHM=2.8 mm. The 
cell was filled with paraffin oil dyed with iodine. The 
iodine served as a weak absorber of the green light and 
its low concentration of 0.5 mg/1 ensured that the total 
absorption in the cell was ^ 10%. The density of paraffin 
oil decreases with increasing temperature, thus resulting 
in a self-defocusing nonlinear response. The two glass 
slides were tilted at a small angle with respect to the 
beam and modified the phase structure of the beam such 
that it results in tt phase jumps in the beam at the posi- 
tion of the glass edge. The phase modulation also gives 
rise to amplitude modulation at the front face of the cell 
as seen in the inset of Fig. ^ The input and the output 
facets of the cell were imaged onto two CCD cameras 
by large numerical aperture lenses. The beam evolution 
along the propagation direction was traced by immers- 
ing a translatable mirror into the liquid and visualizing 
the beam profile at different distances 29]. Because of 
the finite size of this mirror, however, the first 19 mm of 
beam propagation could not be directly accessed. 

The input beam phase profile was monitored by inter- 
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FIG. 4: (color online) Measured output distance between 
two dark solitons as a function of their initial separation. 
Dots - nonlinear regime (3.5 W); Squares - linear propaga- 
tion (10 mW). For small initial separations the output dark 
notches are practically not detectable in the linear regime. 
Inset - interferogram reflecting the initial phase and intensity 
profile of solitons separated by 59 /xm. 



ference with a reference plane- wave [Fig. ^top inset)]. 
Both phase jumps were set to tt within the accuracy 
of our interferometric measurement (~10%). This ac- 
curacy was increased by monitoring the beam profile at 
the far field and setting the intensity in the dark notch 
to zero. By varying the transverse position of the slides 
we could change the initial separation between the gen- 
erated dark notches. A typical intensity profile at a 
separation of 59 /im is shown in [Fig. ^bottom inset)]. 
The width of each notch is a ~ 18 /im, which corre- 
sponds to ~10 diffraction lengths of propagation in the 
nonlinear medium. The individual propagation in the 
self-defocusing medium results in formation of a "black" 
soliton of zero transverse velocity. The experimental con- 
ditions for dark soliton formation were determined by 
attaining saturation of the, so-called soliton constant, 
/oa^, with input beam power 30]. In our experiments 
this regime was reached at approximately 3 W of input 
power. All our measurements were performed at slightly 
higher power levels (3.5 W) to assure soliton regime in 
the presence of interaction and absorption. 

When placed sufficiently close, both dark solitons (with 
parallel initial trajectories) interact during propagation. 
This influences their output separation, which is com- 
pared to the separation at low power. In Fig. ^ we de- 
pict the spacing between the dark solitons at the exit 
of the cell as a function of their initial separation. The 
dots represent nonlinear regime (3.5W), while squares 
correspond to linear propagation (10 mW). The measure- 
ments quite faithfully reproduced the theoretical predic- 
tions. For large initial spacings (> ISO jam) the solitons 
are weakly interacting, but their separation is affected by 
the background beam nonlinear broadening, resulting in 
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FIG. 5: (color online) (a) Measured output separation of 
interacting nonlocal dark solitons in thermal medium as a 
function of the propagation distance for input separation of 
117 /im - closed circles; 101 /xm - open circles; 91 /xm - open 
triangles; 59 /xm - closed triangles, (b) Measured trajectories 
of nonlocal solitons for initial separation of llT/xm. 

increasing separation, which becomes larger than in the 
linear regime. For initial separations of 20-180/im the 
final spacing between the solitons is less than in the lin- 
ear regime. This unambiguously indicates the presence 
of attractive forces counteracting the natural repulsion of 
solitons. Moreover, for input spacings of 115-200 /im the 
final separation is actually smaller than the initial one. It 
should be noted that the separations recorded in the non- 
linear regime are detected with high accuracy (~ 1/im) 
due to the localization of the dark notches. In the lin- 
ear regime the determination of the central position of 
the dark notch is a subject of broadening and reshaping. 
The comparison of these two cases, however, gives us the 
correct information for the range of strong attraction. 

A clearer picture of soliton interaction can be obtained 
by following their trajectories inside the medium. This 
is illustrated in Fig. [HI where we show the measured sep- 
aration along the cell for initial spacings of 59, 91, 101, 
and 117 jam. For the separation of 59 jam the repulsive 
force is strong and can not be compensated by the nonlo- 
cality induced attraction. Therefore, for this separation 
the two dark solitons repel. The situation changes dras- 
tically for larger initial separations. At a certain distance 
the solitons actually come closer than their initial sepa- 
ration. This behavior cannot be directly observed for the 
initial separations of 91 and 101 /xm since the decrease 
in soliton separation occurs at the initial non-accessible 
part of propagation. However, it is particularly visible 
for 117 /im separation. The non- monotonous character 
of the trajectories is a direct proof of the interplay of re- 
pulsive and nonlocality-mediated attractive forces acting 
between the solitons. When the attraction dominates, 
solitons decrease their mutual separation until the repul- 
sion prevails forcing them to move apart. The contour 
plot in Fig. IStb) shows the experimentally obtained tra- 
jectories in this regime (separation of 117 jam). Dashed 



lines indicate location of the intensity minima. The in- 
ward bending of these trajectories is clearly visible. 

In conclusion, we have shown, what we believe is the 
first experimental demonstration of the nonlocality me- 
diated attraction of dark spatial solitons. Our experi- 
mental observations are in good agreement with direct 
numerical simulations. We believe, that our results may 
be applicable to other physical systems exhibiting nonlo- 
cal nonlinear response. 
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